size <- 12. Later this will be the number of rows of the matrix.x <- rnorm( size ).x1 by adding (on average 10 times smaller) noise to x: x1 <- x + rnorm( size )/10.x and x1 should be close to 1.0: check this with function cor.x2 and x3 by adding (other) noise to x.size <- 12
x <- rnorm( size )
x1 <- x + rnorm( size )/10
cor( x, x1 )
[1] 0.9881742
x2 <- x + rnorm( size )/10
x3 <- x + rnorm( size )/10
x1, x2 and x3 column-wise into a matrix using m <- cbind( x1, x2, x3 ).m.m.heatmap( m, Colv = NA, Rowv = NA, scale = "none" ).m <- cbind( x1, x2, x3 )
class( m )
[1] "matrix" "array"
head( m )
x1 x2 x3
[1,] 0.1162679 0.11441132 0.05877585
[2,] -0.1011297 -0.06090018 0.09986234
[3,] -0.8266521 -0.72764170 -0.65511617
[4,] 0.6729362 1.02150546 0.88394058
[5,] 0.2275016 -0.16038791 0.01993610
[6,] -0.1364853 -0.03543970 0.01109748
heatmap( m, Colv = NA, Rowv = NA, scale = "none" ) # high is dark red, low is yellow
# x1, x2, x3 follow similar color pattern, they should be correlated
y1…y4 (but not correlated with x), of the same length size.m from columns x1…x3,y1…y4 in some random order.y <- rnorm( size )
y1 <- y + rnorm( size )/10
y2 <- y + rnorm( size )/10
y3 <- y + rnorm( size )/10
y4 <- y + rnorm( size )/10
m <- cbind( y4, y3, x2, y1, x1, x3, y2 )
heatmap( m, Colv = NA, Rowv = NA, scale = "none" ) # high is dark red, low is yellow
cc <- cor( m ) to build the matrix of correlation coefficients between columns of m.round( cc, 3 ) to show this matrix with 3 digits precision.cc <- cor( m )
round( cc, 3 ) #
y4 y3 x2 y1 x1 x3 y2
y4 1.000 0.994 0.345 0.994 0.365 0.340 0.997
y3 0.994 1.000 0.311 0.995 0.335 0.313 0.993
x2 0.345 0.311 1.000 0.272 0.975 0.994 0.300
y1 0.994 0.995 0.272 1.000 0.301 0.270 0.995
x1 0.365 0.335 0.975 0.301 1.000 0.982 0.321
x3 0.340 0.313 0.994 0.270 0.982 1.000 0.294
y2 0.997 0.993 0.300 0.995 0.321 0.294 1.000
heatmap( cc, symm = TRUE, scale = "none" )
# E.g. value for (row: x1, col: y1) is the corerlation of vectors x1, y1.
# Values of 1.0 are on the diagonal: e.g. x1 is perfectly correlated with x1.
# Correlations between x, x vectors are close to 1.0.
# Correlations between y, y vectors are close to 1.0.
# Correlations between x, y vectors are close to 0.0.
Copyright © 2022 Biomedical Data Sciences (BDS) | LUMC